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Abstract. We consider a three dimensional family of filters based on broken 
power-law spectra to search for gravitational wave stochastic backgrounds in 
the data from Earth-based laser interferometers. We show that such templates 
produce the necessary fitting factor for a wide class of cosmological backgrounds 
and astrophysical foregrounds and that the total number of filters required to 
search for those signals in the data from first generation laser interferometers 
operating at the design sensitivity is fairly small. 



1. Introduction 

Gravitational wave (GW) stochastic backgrounds represent one of the classes of signals 
that could be detected by Earth-based laser interferometers. Stringent upper-limits 
and/or detection will provide new insights into the evolution of the Universe at early 
cosmic times and high energy and, possibly, information about large populations 
of faint unresolved astrophysical sources at low-to-moderate redshift (see e.g. [1] 
and references therein for a recent review). The energy and spectral content of a 
stochastic background of gravitational waves are encoded in its spectrum f2(/) = 
(1/pc) (dp(f)/dlnf), where p(f) is the energy density in GWs and p c is the total 
energy density to close the Universe today. The current analysis of the data from 
Earth-based laser interferometers is restricted to stochastic signals whose spectrum is 
constant over the observational window, i.e. £l(f) — const. However, some theoretical 
models suggest that a wide class of stochastic backgrounds are actually described 
by a frequency dependent spectrum, peaked at some characteristic frequency (that 
depends on the model parameters), with 0(/) following roughly a "broken power- 
law" behaviour; this is true both for cosmological backgrounds and astrophysical 
foregrounds (see for instance [2, 3]). In the near future it is therefore important 
to generalise the analysis of the data for this class of signals. 

In this paper we evaluate whether the use of templates corresponding to a broken 
power-law spectrum f2(/) is suitable to search for stochastic backgrounds characterised 
by a fairly general frequency dependent spectrum. In particular, for this class of 
templates we investigate the fitting factor and estimate the related computational 
costs. 
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2. Filter family and fitting factor 

We consider a family of niters corresponding to a class of broken power-law spectra 
described by 3 parameters: a "knee" frequency f s and two slopes a and b. We choose 
the following functional form for the spectrum 

where fio is simply a normalisation factor and does not represent a search parameter 
and the functions A± are step functions. 

The real signal will most unlikely match exactly the functional form (1), and we 
need to explore the Fitting Factor (FF) [4] of the template family (1). Lacking any solid 
prediction about the expected signals, we have computed the FF of broken power-law 
templates for a broad class of spectra We have considered a Gaussian-shaped 

spectrum, described by the two free parameters /g and A/g, a Lorenzian-shaped 
spectrum, described by the two free parameters /l and A/l, and a broken power- 
law spectrum modulated by an oscillatory term, whose free parameters are, beside 
(a ,b, f s ), the amplitude A 

«g(/;A) 
^osc(/; A) 

The above functional forms Qg(/), and O osc (/) capture the behaviour of the 

spectra that characterise a wide class of cosmological and astrophysical models (see 
e.g. [2, 3]). 

We have computed the fitting factor of the templates £1bpl(/), Eq. (1), for Qg(/)> 
and O osc (/). The results are summarised in Figure 1. It is clear that by suitably 
choosing the parameter space for the slopes a and b one can detect essentially any 
spectrum characterised by a Gaussian/Lorcnzian shape. Even if the real signal shows 
fairly strong and rapid oscillations superimposed to the general power-law behaviour, 
the templates ^bpl(/) can still have a fitting factor greater than 0.97%, except in 
extreme cases. The main outcome of this analysis is therefore that power-law templates 
are indeed suitable for searching for fairly general frequency dependent stochastic 
backgrounds. 

3. Number of filters 

We can now address the computational requirements to set up a search using the 
templates Obpl(/)- In order to compute the total number of filters we have extended 
the geometrical formalism developed in [6] to the case of stochastic backgrounds. The 
basic idea is the following: the signal is a vector in the vector space of the data and 
the A^-parameter family of templates traces out an iV-dimensional template manifold. 
The parameters themselves are coordinates on this manifold, and one can introduce a 
N x N non-flat metric jjk which is related to the fractional loss in the signal-to-noisc 
ratio when there is a mismatch of parameters between the signal and the filter. The 
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Figure 1. The fitting factor of "broken power-law" template families. Left: The 
plot shows, as a function of |a max | = |fo m ax|, the minimum value of A/g (solid line 
top) and A/l (dottcd-linc bottom) corresponding to FF > 0.97 for backgrounds 
characterised by a spectrum (2) and (3), respectively. 

Right: The plot shows, as a function of / OS c, the maximum amplitude A OS c for 
a signal (4) for which FF > 0.97 over a frequency band 40Hz-lkHz. The two 
curves correspond to a knee frequency f s in the signal corresponding to 100 Hz 
(top) and 1kHz (bottom), respectively. The slopes for both f!BPL(/) ar >d f^osc 
are a = 2 and b = —2. 



spacing of the grid of filters is determined by the fractional loss due to the maximum 
tolerable mismatch fi, which fixes the grid spacing of the filters in the parameter space 
V . For a a hyper-rectangular mesh (which does not represent necessarily the most 
efficient tiling of the parameter space) the number of filters M is given by 

N 

i In 
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AJ det|| 7 , fc ||; (5) 



In particular, by adopting a detection strategy based upon the cross correlation 
statistics [5], for a stochastic signal characterised by the spectrum (1) the metric jjk is 
represented by a 3 x 3 symmetric matrix whose components depend on the following 
inner product (involving the optimal filter Q) and its I s * and 2 nd derivatives with 
respect to the slopes (o, b) and to the "knee" frequency f s (a more detailed derivation 
can be found in [7]) 



Jf mt n P Sl(f) S 2 (f) 

where {f mini f max) are the low-frequency and high-frequency cut-off of the 
interferometer band (chosen to be 40Hz and IkHz respectively), Si, 2 are the noise 
spectral densities of the detectors and 7 is the overlap reduction function. 

We have applied this formalism to stochastic backgrounds whose spectrum is given 
by (I) and we have computed N for N = 3 and [i = 0.03 (a 3% maximum tolerable 
mismatch) for the case of cross-correlations between the Livingston and Hanford LIGO 
interferometers operating at the designed sensitivity over 10 7 sec of integration time. 
Figure 2 summarises the results. In particular, we have considered both the more likely 
case of broken power law spectra with a peak in the sensitive region of the detector 
frequency band - e.g. due to a phase transition in the Early Universe - (dash-dot line) 
and the case when also spectra with a "valley" are taken into account (solid line). It 
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Figure 2. The number of niters for a search based on "broken power-law" 
templates. Left: The plot shows the total number of filters, assuming to know 
a priori the frequency f 3 , required to carry out the analysis over 10 7 sec by 
cross-correlating the data streams from the two LIGO interferometers at designed 
sensitivity, as a function of f s . The range of the slopes is —10 < a, b < 10. Notice 
that, as expected, the largest number of filters corresponds to f s right at the 
heart of the sensitivity window. Right: The plot shows the total number of filters 
for the general case of a 3— dimensional parameter search (both on the spectral 
slopes and the knee frequency) . We show the total number of filters as a function 
of the maximum value (a max ) of the slope parameter a, assuming a maximum 
signal-to-noise ratio loss of 3%, for two different choices of the parameter space: 
(i) a m in = bmin = -10, 6 max = a m ax (solid line), and (ii) a min = 6 max = 0, 
frmin — — a max (dashed-dotted line). 

is clear that for a 3-dimensional search over the parameter space 40 Hz < f s < 1 kHz 
and —10 < a, b < 10 - which essentially covers the relevant parameter space (see 
previous section)- the total number of filters that are required for a 3% mismatch is 
^ 500. This represents a modest number of filters, and the search does not pose any 
significant challenge as far as computational time is concerned. 

4. Conclusions 

We have considered a family of templates characterised by a broken power-law 
spectrum to search for fairly general classes of frequency-dependent gravitational wave 
stochastic backgrounds. We have shown that templates corresponding to broken-power 
law spectra have the necessary fitting-factor (i.e. greater than 97%) for a wide class 
of models and that the total number of filters required to carry out a search using the 
data from first generation laser interferometer is sufficiently small not to produce a 
large computational burden. 
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